Jaynes-Cummings photonic superlattices 



Stefano Longhi 

Dipartimento di Fisica, Politecmco di Milano, Piazza L. da Vinci 32, 1-20133 Milano, Italy 



o 

(N 

> 
O 

in 



^ ■ 

Oh- 



> 



X 
S3 



Compiled November 16, 2011 

A classical realization of the Jaynes-Cummings (JC) model, describing the interaction of a two-level atom 
writh a quantized cavity mode, is proposed based on light transport in engineered waveguide superlattices. The 
optical setting enables to visualize in Fock space dynamical regimes not yet accessible in quantum systems, 
providing new physical insights into the deep strong coupling regime of the JC model. In particular, bouncing 
of photon number wave packets in Hilbert space and revivals of populations are explained as generalized Bloch 
oscillations in an inhomogeneous tight-binding lattice. © 2011 Optical Society of America 
OCIS codes: 230.7370, 020.5580, 000.1600 



Transport of discretized light in photonic lattices has 
received a great attention in recent years, suggesting 
the possibility to realize novel functionalized optical ma- 
terials [1, 2]. Photonic lattices also provide a beauti- 
ful laboratory system to realize classical analogues of a 
wealth of quantum phenomena ranging from solid-state 
physics [1-10] to relativistic quantum mechanics [11-15]. 
In this Letter it is shown that photonic lattices can pro- 
vide classical realizations of cavity quantum electrody- 
namics (QED) systems as well, offering the possibility to 
access and visualize dynamical regimes not yet accessi- 
ble in QED systems. Here a photonic realization of the 
renowned Jaynes-Cummings (JC) model, describing the 
interaction of a two-level atom with a quantized mode of 
a cavity [16, 17], is proposed. The JC model has seen a 
renewed interest since the advent of circuit QED, where 
on-chip superconducting qubits and oscillators play the 
roles of two-level atoms and cavities. In such systems, 
the possibility to explore novel dynamical regimes has 
been suggested, and the first experimental observations 
of breakdown of the rotating wave approximation (RWA) 
have been reported [18-21]. Recently, the deep strong 
coupling (DSC) regime of the JC model has been theo- 
retically investigated [20] , and bouncing of photon num- 
ber wavepackets in Hilbert space, leading to collapse and 
revivals of the initial populations, have been predicted. 
The DSC regime, however, is still far from being acces- 
sible in semiconductor or superconducting circuit QED 
systems. The classical realization of the JC model pro- 
posed in this Letter provides a testbed to mimic in the 
lab all the dynamical regimes of the JC Hamiltonian 
directly in Hilbert space. It also gives a new physical 
view onto the revival dynamics of the DSC: it can be ex- 
plained as a kind of generalized Bloch oscillations of the 
wave function in an inhomogeneous tight-binding lattice. 

The JC Hamiltonian describing a qubit interacting 
with a bosonic mode [Fig. 1(a)] reads [17] 

H — + tiuja^a + fig{(j+ + (J_)(d -I- a^), (1) 

where a and are the annihilation and creation op- 
erators of the quantized oscillator with frequency w, 



CTz = |e)(el - \g){gl &+ = le)(5l and ct- = \g){e\ are 
Pauli operators associated to the qubit with ground state 
\g), excited state ]e), and transition frequency wq, and g 
is the coupling strength. The usual JC model under the 
RWA is obtained in the weak-coupling and near-resonant 
limits (]w — wg] u!q, u! ^ g), and corresponds to ne- 
glect the anti-resonant interaction terms a^a) and a_a 
in the Hamiltonian (1). On the other hand, the DSC 
limit corresponds to a coupling g of the order or larger 
that the oscillator frequency w, with no particular re- 
lation between uj and wq [20]. The main idea underly- 
ing the proposed photonic realization of the JC model 
is that the temporal evolution of the quantum system 
in Hilbert space (spanned by the states \g)\n) and |e)ln) 
describing n quanta of the field with the atom in the 
ground or in the excited state) can be mapped into 
the spatial propagation of light waves in a semi-infinite 
curved binary array of waveguides with non-uniform cou- 
pling rates. In fact, let us expand the vector state of 
the system as \il}{t)) = Yln=o\-°'ri{t)\e)\n) + hn{t)\g)\n)] 
and let us introduce the amplitudes c„(t) and fn(t), 
defined by c„ (t) = a„ , /„ (t) = bn for n even, and 
c„(t) = bn , fn(t) = a„ for n odd. Then, the tempo- 
ral evolution of the amplitudes c„ and /„, as obtained 
by the Schrodinger equation ifidt\ij{t)) = H\ip{t)), read 

l—r- = K„C„+i -I- K„_iC„_i + (-1) —Cn+nUJCn (2) 



dt 

dfn 



1^0 



= K7i/ri-Hl + tn-i/n^i - (-1)"— /„ + nw/„ (3) 



(n = 0, 1, 2, 3, ...), where we have set k„ ~ gy/n+l. In 
the optical context, Eqs.(2) and (3) describe light trans- 
port in two uncoupled semi-infinite binary photonic lat- 
tices in the tight-binding approximation with a super- 
imposed transverse index gradient and with non-uniform 
coupling constant k„ between adjacent waveguides n and 
(n-|-l) (see, for instance, [7]). In the optical analogue, the 
temporal variable t plays the role of the spatial propaga- 
tion distance, the transition frequency uq of the qubit 
corresponds to the propagation constant mismatch of 
alternating waveguides of the superlattice, and the fre- 
quency Lu of the quantized mode to the transverse index 
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Fig. 1. (Color online) (a) Schematic of a qubit interact- 
ing with a single bosonic mode, (b) Photonic realization 
of the JC model in Hilbert space based on a semi-infinite 
curved binary waveguide array. 

gradient. The superlattices (2) and (3) are decoupled and 
basically correspond to the two uncoupled parity chains 
of the Hilbert space discussed in Ref. [20] . Since the su- 
perlattices (2) and (3) are obtained each other by revers- 
ing the sign of wq , the analysis can be limited to consider 
the dynamics of either one of the two lattices. Their opti- 
cal realization can be obtained using a semi-infinite chain 
of circularly-curved waveguides with alternating propa- 
gation constants and with engineered distances d„ in the 
geometrical setting depicted in Fig. 1(b) [22]. The index 
gradient lu is given by w = 2TTnsa/ (i?A), where is the 
refractive index of the substrate, A is the wavelength of 
light, R the bending radius of curvature, and a is their 
horizontal spacing [22]. It should be noted that in the 
limiting case w = wq = the lattice reduces to the one 
recently proposed in Ref. [23] as a classical analogs to 
quantum coherent and displaced Fock states. 

The usual RWA of the JC model corresponds to 
strongly detuned waveguides, g/tuo — > 0, with the reso- 
nance condition a; ~ wq, which decouple light dynam- 
ics of the array in pairs of waveguides (i.e. in a se- 
quence of optical directional couplers). In fact, if we 
consider for example the lattice (2), after setting c„ — 
On eyi\)[—i{—l)'^ujQt/2 — inojt], in the RWA approximation 
from Eq.(2) a set of decoupled equations are obtained 
for the slowly-varying amplitudes On and 9n+i describing 
the occupation amplitudes of states |e)|n) and |5)|n-|-l), 
respectively 



l—j^ = Knexp{-lAt)On+l 



n+1 



dt 



where A = lj — ujq and the integer n is taken to be 
even. Hence, Rabi-like oscillations of the optical power 
between the two waveguides n and {n + 1) with spatial 
frequency fin = \/g^(n + 1) -f (A/2)^ is obtained, which 
is the analog of the QED Rabi frequency [17]. 

A more exotic dynamics is observed in the DSC regime 
of the JC model [20] . In the limit of degenerate qubit lev- 
els, i.e. for wo/w —> 0, the energies and eigenstates of the 
JC Hamiltonian can be determined analytically, and the 
dynamics of populations and photon numbers is strictly 
periodic with period T = 27r/w [2]. As an example, in 




Fig. 2. (Color online) (a) Behavior of populations Pg 
(solid curve) and (dashed curve) of ground (|.g)) and 
excited (|e)) qubit states versus normalized time ujt for 
the JC Hamiltonian (1) in the DSC regime for initial 
state |^/'(0)) = |5)|0). Parameter values are g/ui = 2 
and ujq/uj = 0. The corresponding behavior of pho- 
ton statistics P{n,t) = |/?i(i)|^ and revival probability 
Prev = l(^'(OIV'(0))|^ are shown in (b) and (c). 



Fig. 2 it is shown the dynamical evolution of populations 
in the two qubit levels, Pg{t) and Pe{t) = 1 — Pg(i), 
and of the photon number distribution P(n, t) as ob- 
tained from numerical analysis of the JC Hamiltonian 
(1) when the system is initially prepared in the pure 
state \il}{t = 0)) = |5)|0). In the figure, the evolution 
of the return (revival) probability to the initial state, 
P^g^(t) = |(V'(i)|V'(0))|^, is also depicted. The behav- 
ior of Prev{t) clearly shows the existence of collapse 
and (exact) revivals, which have been explained in Ref. 
[20] as due to the periodic bouncing of photon number 
wavepackets along one of the two parity chains in the 
Hilbert space, as shown in Fig. 2(b). Such a behavior is 
deeply distinct from the usual collapse-revival dynamics 
of the JC model in the RWA, which requires large initial 
coherent states [17]. In our optical realization of the JC 
model, the evolution of the light beam intensity along the 
array of Fig. 1(b) when the boundary (n = 0) waveguide 
is excited at the t = input plane exactly reproduces 
the evolution of the photon number distribution shown 
in Fig. 2(b). In fact, the populations Pg{t) and Pe{t), the 
photon number distribution P{n, t) and revival probabil- 
ity Prev {t) are related to the occupation amplitudes fnit) 
of light beam in the waveguides of the lattice according 

to Pgit) - Er=ol/2n(OP, Peit) = E^^O I (0 P , 

Pin,t) = |/„(^)|^ and Previt) = |/o (<)!'• I* should be 
noted that the proposed optical analog simulator gives 
a new physical view of the revival dynamics of the JC 
model in the DSC regime not noticed in previous work 
[20] . Indeed, the periodic bouncing of the photon wave 
packet in Hilbert space shown in Fig. 2(b) is the signature 
of the self-imaging property of the waveguide lattice, and 
the observed oscillatory beam path can be ultimately ex- 
plained in terms of a kind of generalized Block oscilla- 
tions in an inhomogeneous tight-binding lattice under a 
dc bias. In fact, in the limit wq/w — >■ the lattice de- 
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scribed by Eq.(3) belongs to a class of exactly-solvable 
lattice models which shows an equally-spaced Wannier- 
Stark ladder spectrum with energies Ei = lui — g'^ /oj 
{I = 0,1,2,3,...), as shown in Ref. [24]. Note that, as 
compared to most common optical Bloch oscillations 
previously investigated in homogeneous lattices [1,4,5], 
in our case the waveguide lattice turns out to be trun- 
cated and inhomogeneous. However, the main oscilla- 
tory dynamics of the optical beam can be still explained 
from a semiclassical analysis of the inhomogeneous lat- 
tice Hamiltonian as a result of beam acceleration induced 
by the dc bias field [the last term on the right hand 
side of Eq.(3)] and Bragg scattering in the lattice. More- 
over, since for luq/uj ^ the JC model is realized by a 
single band tight-binding lattice, Zencr tunneling [6, 7] 
or other phenomena, like interband Rabi oscillations of 
Bloch modes, do not occur in our case. 

Typical design parameters of the lattice that realizes 
the dynamics of Fig. 2 can be obtained for femtosecond 
laser written waveguides in fused silica [2] following the 
lines detailed in Ref. [22]. For instance, if we assume cir- 
cular waveguides with core diameter p = 5 fim and index 
change An = 0.002 (as in Ref. [22]), for an axis bend- 
ing radius i? = 60 cm, an horizontal waveguide spacing 
a = 6 /xm and a bulk refractive index = 1.45, at 
the wavelength A ~ 633 nm the period of Bloch oscil- 
lations (i.e. of the revival) shown in Fig. 2(b) turns out 
to be T = 27r/a; ~ 4.37 cm. The coupling constants arc 
given by k„ = g\/n + 1 with 9 = 2^^ 0.288 mm^^. 
The distances do, di, d2, ... between adjacent waveg- 
uides [see Fig. 1(b)] that realize such non- uniform cou- 
pling constants are calculated by the approximate rela- 
tion [22] dn ~ (l/7)ln(A/K„) with A ~ 24.6 mm"! and 
7 ~ 0.466 /im^^, i.e. do — 9.54 ^m, di ~ 8.80 ^m, 
d2 — 8.37 //m, etc. Note that a rather limited number of 
waveguides, less than ^ 25 [see Fig. 2(b)], is enough to 
visualize the dynamics in the Hilbert space. Finally, it 
should be noted that a nonvanishing value of wq/w makes 
the revivals only approximate, as shown in Fig. 3. In this 
case, the lattice is a binary one and the Wannier Stark 
states are no more equally space. However, even though 
the wave packet revival is incomplete, photon bouncing 
is clearly observable over a few cycles. 



(a) (b) 




Fig. 3. (Color online) Same as Fig. 2, but for g/uj ~ 2 
and wo/w = 0.3. 



In conclusion, an optical realization of the JC model 
of QED, based on light transport in waveguide super- 
lattices, has been proposed. Our optical setting enables 
to simulate the physics of atom-field interaction in the 
deep strong coupling regime [20] , which is currently not 
accessible in cavity or circuit QED systems. 
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